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Abstract 

In this talk, I review the Color Glass Condensate theory of gluon saturation, 
and its application to the early stages of heavy ion collisions. 
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1. Color Glass Condensate 

In heavy ion collisions, the vast majority of the produced particles have a 
transverse momentum of at most a couple of GeV's. This implies that they 
are produced from incoming partons that carry a very small fraction x of 
the longitudinal momentum of the projectiles. Furthermore, this momentum 
fraction decreases with the center of mass energy of the collision. At RHIC 
energy, the typical value of x is of the order of 10~ 2 for the bulk of particle 
production, and of order x ~ 4.10 -4 at the LHC. However, at such small 
values of x, the gluon density in a proton or nucleus becomes large which 
leads to a new phenomenon known as gluon saturation, due to non-linear 
recombinations among the gluons. Gluon saturation is characterized by an 
x-dependent momentum scale, the saturation momentum Q s (x), that delim- 
itates the domain in momentum where saturation takes place. This domain 
is represented in the figure [TJ 

Having identified this region where nonlinear gluon interactions become 
important, one faces now the question of how to compute physical observ- 
ables reliably in this regime where one collides projectiles made of a very 
large number of constituents. The main novelty in the saturated regime, 
compared to the dilute one, is that the typical particle production processes 
involve multiple gluon interactions, as shown in the figure |2j The standard 
perturbative techniques based on Feynman diagrams are well suited to han- 
dle the dilute situation, but they become impractical in the saturated regime 
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Figure 2: Left: Gluon production in the dilute regime. Right: Saturated regime. 



because an infinite number of graphs would contribute at each order in g 2 . In 
addition, in order to compute observables in this dense regime, we also need 
to know the probability to find these multigluon states in the wavefunctions 
of the incoming nuclei. 

The Color Glass Condensate is an effective description of QCD in the 
saturation regime. In this effective theory, one divides the degrees of freedom 
(see the figure [3]) into fast partons with longitudinal momentum k + > A + , 
and slow partons with k + < A + that have a significant dynamical evolution 
over the time-scales of the reaction of interest jlj]. Thanks to Lorentz time 
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Figure 3: Separation of the degrees of freedom in the CGC. 

dilation, the fast partons are essentially frozen during the collision process, 
and it is sufficient to know what color current they carry. Therefore, their 
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description in the CGC is simplified into a distribution of static color currents 
jfj. = 5 At+ p(x~, x±) (this is for a projectile moving in the +z direction). 
On the contrary, this approximation does not work for the slow degrees of 
freedom, and therefore we must keep describing them as conventional gauge 
fields A 11 . Due to the separation in longitudinal momentum between these 
two types of degrees of freedom, their coupling is eikonal, A^J^. Thus, the 
effective Lagrangean of the CGC is: 

C= —txF^F^ + (Jf + J£) A,. 

gluon dynamics fast partons 

This description would not be complete without a way to specify the color 
charge density p(x~,x±). This distribution depends on the precise configu- 
ration of the fast partons (e.g. their localization in the transverse plane at 
the time of the collision), which is not known. Thus, it should be considered 
as a stochastic quantity, with a probability distribution Wa+[p]- The distri- 
bution Wa+[p] plays in the CGC the same role as partons distributions in 
conventional low density perturbative QCD. 

The CGC effective theory introduces a cutoff A + that separates the fast 
and the slow degrees of freedom. However, this cutoff is unphysical and one 
should ensure that physical observables do not depend upon it. It turns out 
that for this to be true, the probability distribution W\+ [p] must evolve with 
A + according to the JIMWLK equation ^ , 

H W A + , H = - / -=—. — -r](x ± ,y ± ) 



<91n(A+) A ' 2 J 5a(y ± ) ,x ^^^'d^x^) 1 

x±,y± 

where — d\ a(x±) = p(l/A + , x±). Effectively, the JIMWLK equation resums 
to all orders the powers of a s ln(A + ) (or, equivalently, powers of a s \n(l/x)), 
and its kernel r)(x±, y ± ) includes all corrections in the color source p. When 
approximated to low density by expanding the kernel to lowest order p, one 
recovers the BFKL equation that describes the evolution to low x of non- 
integrated gluon distribution of a dilute hadron. 

2. Just before the collision: factorization 

A high energy heavy ion collision can be divided in several successive 
stages, as shown in the figure [TJ The CGC applies to the description of the 
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wavef unctions of the incoming projectiles, to the collision itself, and to a brief 
time after the collision, of the order of r ~ Q' 1 . Since the subsequent stages 
are usually described by some form of relativistic fluid dynamics, one of the 
CGC goals is to provide the value of the energy momentum tensor at some 
initial time To- But more fundamentally, since the CGC provides a QCD- 
based description of the early stages of heavy ion collisions, it is also the 
framework to be used in order to justify the applicability of hydrodynamics. 

In the CGC, the energy-momentum tensor admits an expansion in powers 
of the strong coupling g 2 . However, because the color sources are strong in 
the saturated regime, this expansion starts by a term in g~ 2 , 

c + cig 2 + c 2 g A H . 

The coefficients £0,1,2,- in this expansion are themselves infinite series in gp 
(which is parametrically of order one in the saturated regime), and therefore 
each order corresponds to an infinite series of Feynman diagrams. 

At leading order, one can prove that the sum of the corresponding set of 
diagrams can be expressed in terms of classical solutions of the Yang-Mills 
equations 0, 

LO 4 1 J / t-»— 00 

Yang— Mills equation 

These equations are non-linear and in general one cannot find analytic solu- 
tions. However, they have been solved numerically to obtain e.g. the initial 
energy density released in a collision Q. 

At next to leading order, the energy-momentum tensor starts being sensi- 
tive to the slow gluon fields of the CGC effective theory, via a loop correction. 
It is convenient to consider only one small slice of these slow modes, located 
just below the original cutoff Aq~ of the CGC, as shown in the figure HI The 
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Figure 4: Slice of slow modes at NLO. 
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contribution of this slice can be calculated at leading log accuracy |5| 




KS) %+in (l) 



where ?ii t 2 are the JIMWLK Hamiltonians of the two projectiles. Thanks 
to their form that does not mix the two projectiles, these logarithms appear 
to be intrinsic properties of the wavefunctions of the two nuclei. Therefore, 
they can be absorbed by defining a new CGC effective theory that has its 
cutoff at the lower values At, 



provided the distribution of sources in the new effective theory is defined as 
prescribed by the JIMWLK equation. This process can be iterated, until the 
cutoffs become smaller than the physically relevant scales. The outcome of 
this procedure is the following formula for the energy-momentum tensor, 



which resums all these logarithms. Note that at this level of accuracy, the 
rapidity dependence of the left hand side comes entirely from the JIMWLK 
evolution of the distributions W of the two projectiles. Besides its useful- 
ness in heavy ion collisions, this factorization is important because it estab- 
lishes a bridge between nucleus-nucleus collisions and other reactions such as 
electron-nucleus collisions, where distributions W with the same JIMWLK 
evolution also appear. Let us also stress that this factorization also applies to 
any sufficiently inclusive observables. In particular, the correlation between 
the components of the energy-momentum tensor measured at different points 
in space [^j reads 



One sees that at leading log accuracy, all these correlations come from the 
Ws (since the integrand is a product of n independent factors). This formula 
predicts a correlation in rapidity over a range of order Ar/ ~ a J 1 , since this 
is the rapidity variation necessary to produce a change in the Ws. 




(T^(r, V ,x ± )) Lhog = J [D Pl Dp 2 ] W 1 [p 1 ] W 2 [p 2 ] T£(t,x ± ) 
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Figure 5: Left: field components evaluated by solving numerically the Yang-Mills equations 
(from [6|). Right: longitudinal color flux tubes. 

3. Just after the collision: Glasma fields 

Immediately after the collision, the chromo- E and B fields have only 
longitudinal components [6], forming flux tubes along the collision axis (see 
the figure |5]). This configuration of color fields has been named the glasma. 
The typical transverse size of a flux tube is of order Qj 1 , and the color fields 
are correlated over a~ l units of rapidity in the longitudinal direction. 

This particular topology of the post-collision color fields has several con- 
sequences, among which a peculiar form of the energy-momentum tensor (see 
section Hj), the fact that the multiplicity distribution is a negative binomial 
0, and the existence of a non-zero topological density FF, possibly at the 
origin of observable CP violating effects jsj. But the most direct and striking 
consequence of these structures, when taken as initial conditions for hydro- 
dynamical expansion, is that they lead to the formation of the so-called ridge 
correlations, a structure in the 2-hadron spectrum which is elongated in Arj 
and narrow in A0 (see the left plot of the figure EJ). By examining the causal 
relation between two particles separated in rapidity (right part of the fig- 
ure [6]), one can see that the process responsible for producing a correlation 
between these particles must have taken place at early times, 

^correlation — ^freeze out ^ 2^ A I _ 

Since the color fields produced at early times in the CGC formalism are 
correlated over rapidity intervals of order a" 1 ^> 1, they provide a natural 
explanation for the rapidity dependence of the ridge [9j. The strength of 
the 2-particle correlation is controlled by (Q S R)~ 2 -the area of one flux tube 
relative to the total transverse area- since the two particles must come from 
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Figure 6: Left: STAR result on 2-hadron correlations. Right: causal relationship between 
two produced particles. 



the same tube to have been produced by the same coherent field (see the 
figured left panel). The azimuthal dependence is produced at a later stage, 
by the radial hydro dynamical flow, that collimates the azimuthal angles of 
the two particles in the direction of the radial velocity (figured right panel). 




Figure 7: Left: particle emitted from distinct tubes are uncorrelated. Right: collimation 
due to the radial flow. 



4. Matching to hydrodynamics 

At times r 3> Q~ l , the standard description of the evolution of the fireball 
is via hydrodynamical expansion. However, a trivial consequence of the fact 
that the chromo- E and B fields are initially parallel to the collision axis 
in the glasma is that the energy-momentum tensor one obtains at leading 
order in g 2 in the CGC is of the form T^(0 + , rj, x) = diag (e, e, e, — e) where 
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e is the energy density. The fact that the longitudinal pressure is negative 
is problematic for a smooth matching to hydrodynamics, because it means 
that the CGC energy-momentum tensor is quite far from the form for which 
quasi-ideal hydrodynamics is applicable (i.e. one where the spatial part of 
the tensor is not too far from being proportional to the identity). 

It has been known for some time that corrections to the leading order 
CGC prediction suffers from secular divergences. Indeed, the solutions of 
classical Yang-Mills equations are usually unstable, as noticed in several 
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This implies that the next to leading order correction to the en- 
ergy momentum tensor calculated in the CGC framework grows with time, 
and eventually becomes larger than the leading order. It is possible to resum 
a subset of these higher order corrections, that at each loop order picks up 
the most singular of these contributions. This resummation amount to su- 
perimposing a fluctuation to the initial condition at r = + for the classical 



glasma field , 



T^ uiamed (r, V ,x ± ) = / [D Pl Dp 2 ] [Da] F[a] W 1 [p 1 ] W 2 [p 2 ] T%[ A + 



initial field 

with a Gaussian distribution F[a] for this fluctuation. The effect of this re- 
summation has not been fully investigated in the context of the CGC due to 
some technical difficulties that preclude a straightforward numerical evalua- 
tion of the previous formula. However, it is possible to devise a scalar toy 
model where the effect of such fluctuations can be studied. To keep things 





Figure 8: Left: evolution of e and p without fluctuations of the initial field. Right: 
evolution with initial fluctuations. 

extremely simple in this toy model, we discard any spatial dependence, so 
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that the classical equation of motion reads 

<j> + -<j>- v 2 x 4> + v'(4>) = o , 

r 

and we take a quartic potential V(<f>) ~ 4 . The term <p/r is due to longitu- 
dinal expansion in the rapidity direction, as would be the case in a collision 
process. The time evolution is started at some some initial proper time 
t = tq, and we average the energy-momentum tensor over Gaussian fluctua- 
tions of the initial field. We first did this toy calculation without fluctuations 
of the initial field (see the left plot of the figure [8]) and there one sees that 
the pressure oscillates with time: there is no single-valued relationship be- 
tween the energy density and the pressure, which means that there is no 
equation of state in this case. In the right plot of the figure [BJ we show the 
effect of the initial field fluctuations: the oscillations of the pressure are now 
damped, and it relaxes towards e/3 - which is the equation of state expected 
for a scale invariant theory in four dimensions. Moreover, one also observes 
that the energy density decreases as r _4//3 , which is the expected behavior in 
hydrodynamics if the longitudinal pressure is e/3. Thus, in this toy model, 
it appears that the fluctuations of the initial field play a crucial role in the 
relaxation of the system towards hydrodynamical behavior. 

Conclusions 

The Color Glass Condensate provides a first principles framework to study 
the early stages of high energy nuclear collisions, that involve gluons in the 
saturated regime. It correctly describes the energy density released in such 
a collision, as well as important features of the correlations observed in the 
final state. In this context, the most pressing and challenging question for 
future work is arguably that of thermalization. 
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